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ABSTRACT 


In the present work, an attempt has been made to use 
mathematical formulations of surfaces in computer-aided- 
design of three-dimensional free-form sculptured surfaces- A 
software package, SURF-3D, has been developed and installed 
using spline, Bezier and B-Spllne representations. Spline 
representations are based on the classical ’lofting' 
techniques. Bezier and B-Spline schemes are easily adaptable 
to interactive design environments because of their Inherent 
properties of ease in definition and manipulation. The 
package helps in designing and displaying surfaces of 
various types. These schemes help in designing axl -symmetr i c 
components also. In addition, alogorithms for viewing 
designed surfaces in parallel or perspective projections, 
for intersections of surfaces with planes and for generating 
hidden surface and hidden line displays have also been 
devised. The package has been tested using several 
Illustrative examples . 
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CHAPTER I 


INTRODUCTION 


1 . 1 Surfaces in Engineering. 

Surfaces and surface represent at i ons play a vital role 
in most engineering design and manufacturing processes. The 
manufacture of aesthetically and aerodynami cal ly designed 
automobile and aircraft bodies, hull design of ships, 
machine parts, design of glassware, clothing, shoes etc. are 
obvious examples. 

For all such surfaces, the present practice of design 
Is based on the use of templates [31. Templates are nothing 
but sheets cut according to a particular contour obtainable 
from a surface when it is cut by a plane at a particular 
position. An ensemble of such templates collectively 
describes the nature of the surface envisaged. 

There are obvious disadvantages to such an approach for 
design and manufacture of surfaces. The overall procedure of 
design, manufacture and modification Is time consuming. 
Visual inspection is a must for such surfaces and therefore 
results in costly prototyping. 

Uith the recent advances in Computational Geometry 
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concepts, it is now possible to mathematically represent 
such surfaces [3,1]. Closed form mathematical equations for 
regular surfaces, such as conics, paraboloids, hyperboloids 
etc. were available even earlier. Use of such concepts has 
given credence to representation of even free-form sculptured 
surf aces . 

Computer Graphics has come in to play a vital role in 
this respect. Nov surfaces can be visualized in their proper 
perspective before manufacture. Analyses of various types 
can be carried out on the mathematical models of such 
surfaces. It has, therefore, proved to be a ’pre-service' 
test . 


Generally surface representations have two different 
uses, 'analytic' and 'synthetic' [8]. 

Analytic uses require a precise ’fit’ of the surface to 
the data points on the surface. Examples of these could be 
turbine blades, aerofoil sections of aircraft, etc. 

Synthetic use of such representations lends itself to 
creation, modification and display of a model of a shape 
interactively, and examination and improvement of the design 
until it is acceptable. The resulting surface may either be 
used to generate a new image or for generating numerically 
controlled machine instructions for actual manufacturing of 
the surface [ 8 ] . 
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1 - 2 Revi «v of developments . 

Various representational schemes have been formulated 
and studied to properly model such surfaces. Earliest among 
these are the classical ’lofting’ techniques of the ship- 
building, automobile and aircraft industries [10]. A Spline 
surface representation was generally used for lofting. 

In the mid 1960s S.A. Coons (1963) of NIT and J.C. 
Ferguson (1964) of Boeing Corporation began work in 
sculptured surfaces [1]- At about this time. General 
Motors developed their DAC-1 system [7]. Several other 
companies, notably Lockheed and NcDonnell Douglas also made 
signi f leant developments [ 7 ] . 

Both Coons and Ferguson used parametric cubics for 
their surface representations. Bezier (1972) evolved a 
successful surface design system UNISURF, which has been 
used in the Renault Co. of France for the design of 
automobile bodies [2]. In the mid 1970s R.F. Riesenfield 
used B-Splines for representation of surfaces for the first 
time [1]. Their emphasis was on the accurate representation 
of surfaces and . on interpolating techniques to meet a 
variety of design criteria. 

The different representational schemes have been 
primarily described in literature on differential geometry, 
computational geometry and computer graphics. Several 
excellent titles such as Faux and Pratt [3], Rogers and 
Adams [10], Mortenson [7] describe the formulations for 
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design of different types of surfaces. An excellent treatise 


on practical use of surface design and manufacturing is by 
Bezier [ 2 ] . 

State-of-the-art techniques in design and display of 
surfaces can be found in Gregory [5] and Barnhill and 
Ri esenf i eld [ 1 ] . 

Viewing and display of surfaces, including shading and 
rendering can be found in Newman and Sproull [8], Foley and 
Van Dam [4], Rogers and Adams [10] and Rogers [9]. Useful 
algorithms for viewing and display can also be obtained from 
Harrington [ 6 ] . 

1 . 3 Surface design methodology. 

In this thesis, the treatment of the mathematical 
techniques applied to modelling of shape is primarily 
concentrated upon the synthetic approach- This approach has 
been adopted because it helps in guiding the designer to go 
through the stages of design and display in a more logical 

i 

and intuitive manner than the analytic approach. Furthermore 
most synthetic methods can be easily applied to analytic 
fitting problems (especially for fitting smooth surfaces), 
whereas analytical methods are not easily amenable to 
problems of free-form and sculptured surfaces. 

The modelling system, therefore can support a class of 
shapes that is matched to the designers' requirements. It is 
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possible to model simple mechanical component shapes, or on 
the other hand, to design as complex a shape as an 
automobile body or an aircraft fuselage. 

The design process in the program is interactive, i.e. 
the designer first visualizes a shape, then applies 
modifications to certain parameters to change the shape, 
until he is satisfied with the design. 

The first step, is therefore, the generation of images 
of shapes that have sufficient realism to communicate to the 
designer the shape information needed to evaluate the 
design. 

The second step is difficult, since the designer has to 
express the required changes in a form which is supported by 
the modelling technique. It has been kept in mind that the 
modifications required by the designer should be expressed 
in such a way that a minimal number of parameters need to be 
changed to achieve the desired result. Here, it is required 
that the designer is sufficiently familiar with the shape 
representation to be able to identify which parameters to 
change and by how much to anticipate very accurately the 
shape changes that will result. 

This last step of selecting parameters to change, 
illustrates the problem of selecting a shape representation: 
if only a small set of parameters is provided, then the 
designers' task is simplified, but the range of shapes that 
can be represented is likely to be limited. If there are 
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many parameters, the designer can achieve more varied and 


more accurate models of his design, but may be overwhelmed 
by the flexibility offered. 


Consequently, keeping these in view, shape 
representations have been chosen which allow a variable 
number of parameters: the designer need use as many or as 
few as necessary to create an acceptable model. 

Mathematically tractable and computationally efficient 
formulation schemes have been adopted in this thesis. These 
representational schemes come close to matching the 
designers’ ’feel’ for design and modification of surfaces. 


The above 
properties of the 


considerations combine to outline 
mathematical representations used 


the 


(i) Control points : 

One of the most common ways of controlling the shape 
of a surface interactively is to specify points through 
which the surface must be fitted or that which control the 
surface shape in a predictable way. 

(ii) Multiple values : 

A three-dimensional surface is usually not a single- 
valued function of two coordinates. The surface may loop 
over itself. 
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(iii) Axes Independence : 

The shape of the surface should be same, no matter what 
coordinate system is used. If the control points are rotated 
then the surface should rotate by the same amount without 
change in its shape. 

(Iv) Global or Local control : 

Uhen a control point Is changed, the shape of the 
surface may be altered only in the region of the control 
point ( local ) , or It may change the shape throughout (global 1 . 

(v) Variation diminishing property : 

Some mathematical representations have a tendency to 
magnify rather than smooth, small irregularities in the 
shape outlined by the control points. Others, possessing a 
variation-diminishing property, always smooth the designers' 
control points. 

(vi) Versatility : 

A shape representation must be flexible enough to 
enable a designer to model a variety of shapes using control 
points as tools. 

Uith these requirements in view, it was decided to 
Include only those representations which fulfil these 
criteria. 

All the representations used in this thesis are based 
upon the parametric, vector-valued function. A point on a 
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surface is represented as a vector -- a function of two 
parameters u and v. (See figure 1 . 1 ) 

T 

r(u,v) = [ x(u,v) y(u,v) z(u,v) } 1.3.1 

As the parameters take on values, usually between zero 
and one, the parametric functions trace out a curve on a 
surface . 

The following three mathematical representations have 
been implemented 

(1) Spline blended surfaces 

(2) Bezier surfaces 

(3) B-Spline surfaces. 

The Spline representation fits a smooth surface 
which passes through the control points. Surfaces, which are 
presently in use may be input using this technique. 
Modification to a control point in this scheme results in 
unpredictable changes, and therefore not suitable for use 
for changes in shape. 

Bezier surfaces, on the other hand only approximate the 
control points. Modifications can be made in an intuitive 
manner. However, a change made to a control point results 
in the alteration of all or major part of the surface, 
thereby making it unsuitable for fine adjustments. 

B-Spline surfaces behave similar to Bezier surfaces. 
The surface approximates the control points closer than a 
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Fig. 1.1 A parametric surface with u & v family of 
curves. 
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Bezier surface. Modifications to a B-Spline surface are 
local in nature. 


Facilities to convert from one surface 
has been provided for different types 
r equireraents . 


type to another 
of modification 


Intersecti one with 
visualized by invoking the 


arbitrary planes may 
intersection commands . 


also 


be 


Display facilities have 


designer to 

vi ew 

the surface 

and perspective 

pro j ect ion 

included in 

the 

viewing op 

manipulat e . 




been provided which enables the 
in any orientation. Parallel 
transformations have been 
erations, which the user can 
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CHAPTER II 


PARAMETRIC CURVES AND SURFACES 


2 . 1 PI f f er ent ial Geometry concepts in curve and surface design. 


A surface defined as a bl-parametric equation is of the 

f orm 


T 

r(u,v) = [ x(u,v) y(u,v) z(u,v) ] 


2 . 1.1 


where u and v are the parameters and c defines the vector to 
the surface, where 


jr = r (u , v ) 


2 . 1.2 


A normal vector n(u, v) to a surface at a point defined by 
r(u,v) is given by 



2 . 1.3 


where n represents the unit normal vector as a result of the 
vector product of the tangent vectors 3r/3u and 9r/5v at the 
point. The order in which the vector product is taken, 
determines the direction of n. (See figure 2 - 1 . 1 ) 
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The equation of t, a plane tangent to a surface at a 
point r is given by 

i ■ frvigj *|5) - o 

where r q Is a vector defining a point on the tangent plane, 
r? is the vector defining the point p on the surface at 
which the plane Is tangent. (See figure 2.1.2) 

2 . 2 Spline curve and surface . 

The equation for a single parametric cubic spline 
segment, in terms of a parametric variable w, is given by 

rfw) - w 0 < w < 1 221 

i-1 i 

T 

where r(w) = [ x(w) y(w) z(w) ] . £_(w) may be considered as 

the position vector of any point on the spline segment. 
[10]. The coefficients are determined by specifying four 
boundary conditions for the spline segment. (See figure 
2 - 2 . 1 ). 

The required boundary conditions for each cubic 
segment consist of two end-points and the tangent vector at 
each end-point. For a single segment between and r^ , 

* 

these conditions are given by 
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2 . 2.2 


ICO) - £l 
£(U - £2 



Therefore the parametric spline segment, making use of 
the boundary conditions, is defined as 

r(w) = r (1 -3w J t2w 3 ) * r (3 w 2 -2w 3 ) ♦ r' (w-2w^ w 3 ) 

1 2 " | 

r - (-w 2 +w 3 ) 2.2.3 

2 


Nov, taking a surface represented by two parameters u and v; 
a spline surface which satisfies the boundary conditions 
along the u edges is given by 

r(u,v) - r(0,v)(l-3u 2 i-2u 3 ) » r(l,v)(3u 2 -2u 3 ) 

^■(0,w).(u-2u 2 +u 3 3 ♦ ^U.vU-u^u 3 ] 2.2.4 

Similarly a spline surface, which satisfies, the 
boundary conditions along w edges is given by 

r(u,v) - r(u, 03.(1 -3v 2 +2v 3 ) + r(u,l).(3v 2 -2v 3 3 

+ f*(u,0).(v-2v 2 -v 3 ) ♦ |£(u,13.(-v 2 + v 3 ) 2.2.5 

Figure 2.2.2 depicts a spline blended surface. 
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2 . 3 Bezier curve and surface . 


The mathematical basis of the Bezier curve is a 
polynomial blending function which interpolates between the 
first and last vertices of a characteristic polygon defining 
the curve. The Bezier polynomial is related to the Bernstein 

polynomial. Thus, the Bezier curve is said to have a 
Bernstein basis. The basis function is given by 


2.3.1 


2.3.2 

with n being the degree of the polynomial and i the 
particular vertex in the ordered set r^ to r\^. In general an 
nth order polynomial is specified by n+1 vertices. The 
points on the curve are then given by 

n 

I(w] = Z L J (w) 0 < w < 1 2.3.3 

i-0 1 n,i 


" n C i .w i .(l-w) n ~* 


" c . -mfe 


The surface equation for a Bezier 
parameters u and v and is algebraically 


surface 
def ined 


involves 

as 


two 
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0 < u,v < 1 


2.3.4 


r(u,v) 


n m 

I 1 l J 
i=OJ =0 J J m > j 


(u).J (v) 

nj 


See figures 2,3.1 and 2.3.2 for Bezier curve and 
surface. 


2 . 4. B-Splin« curvi and surf *c+ ■ 

B~spline curves are similar to Bezier curves in that a 
set of blending functions combines the effects of n+1 
control points rj . The order of these polynomial blending 
functions is specially controlled by a parameter k. The 
b_spline basis functions are recursively defined as 



1 

0 


if tj < w < 
otherwise 


Vl 


2.4.1 


and 


N Cu) 


i,k 


Cw-tjJNjjt-iCw) 

fi+k-l _t i 




2.4.2 


where k controls the order of the ■> resulting basis 
polynomials in the parameter w and thus also controls the 
continuity of the curve. 

The ti are called knot values and they relate the 
parameter w to the r^ control points. 

1 A 



The points on a b spline curve are then given by 



r(w) - 

n 

N (w) 

1=0 i i,k 

2.4.3 

where 





- 0 

if i < k 



- i-k+1 

If k < i < n 

2.4.4 

t i 

* n-k+2 

If i > n 


with 

0 < 

i n+k 

2.4.5 

and 

0 < 

w < n-k+2 

2.4.6 

For 

b-spl ine 

surfaces, extending 

equation C2.4.3) 


take two parameters u and v, 

m n 

r(u,v) - Z 2 1 N Cu).N (v) 2 4 7 

1=0J=0 ij i,k J.J Z - 4 - 7 

where m and n define the control net of points and k and 1 
indicate the orders of the polynomials in the u and v 
directions respectively. See figures 2.4.1 and 2.4.2. 
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2 . 5 Rational parametric curves and surf aces . 

There are certain advantages to be gained by 
defining curves and surfaces in terms of functions more 
general than ordinary polynomials as used in the 
formulations described in the earlier sections . 

Conic sections can only be approximated by the 
polynomial blending functions used in spline, Bezier and b- 
spline formulations. Rational parametric functions allow 
exact representations of such curves and surfaces. 

It has been shown that a parametric polynomial curve of 
degree n can be written in terms. of its basis function B](t) 
as (see equations 2.3.3 and 2.4.3) 
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and wj are scalar weights. Note that if the B^ ( t ) sum to 

one (as is the case with spline and Bezier curve segments ). 

the above equation reduces to simple polynomial form. The 

rational form is therefore more general; the weights may be 

manipulated to gain greater freedom in curve and surface 

definition. The rational form can be re-expressed as a 

polynomial in terms of four-dimensional vectors using 

homogeneous coordinates. On multiplying throughout by w(t) 

equation (2.5.2) and writing g , i) and T for the 

I i 1 

components of r • , the component equations 

w(t)x(t) - 254 

wCt)y(t) * 

w(t)z(t) - ^ i v/ i B i [t) 

when taken together with equation (2.5.3), these equations 
may be summarized as 


get) - ig b (t) 
i=0 1 i 


2.5.5 


in which 


Ei - [ W,*J WjTtj Wj r f Wj ] T 


2.5.6 


and 


2.5.7 


R 


wx wy wz w ] T 


Rational parametric surfaces require a similar 
extension using two parametric variables u and v. Further 
discussion is given in Faux and Pratt [3]. 


2 - 6 Bezier control points o f a surface . 

Starting from a spline-blended surface 'fitted* over a 
m x n grid of control points, it is possible to generate the 
Bezier control points which interpolates the same surface. 

The advantage of using this technique is mainly to help 
in later modification of the surface in a much more 
intuitive and easier method than by changing the positions 
of the surface points themselves. 

Parametric spline - blended surfaces have their 
parametric intervals 0 < u <_ 1 and 0 _< w <_ 1 for each patch 
on the surface, whereas Bezier surfaces have the same 
interval over the entire surface. Therefore 
r eparametr izat ion of the spline-blended surface is required. 
Since the number of points in the u-direction is m+1, an 
obvious choice for the stepsize of the parametric variable u 
is 1/m. Similarly in the v-direction the stepsize is 1/n. 

A Bezier surface is as given in equation (2.3.4), where 
r(u,v) is the resulting surface. Rewriting equation (2.3.4) 


23 


in the matrix form 


[ J(u) ] . [r] . [ J(v) ] = [ r ( u , v ) ] 2.6.1 

where [r] is the unknown matrix of control points for the 

Bezier surface, and [J(u)] and [J(vl] are the square 

matrices of blending function values in the u and v- 

directions respectively, of the Bezier surface. 

[ J(u) ] and [ J(v) ] are positive definite, and hence 

their inverses exist. Premultiplying equation (2.6. I - ) by 
-1 

-1 

[r].[J(v)] = [ J(u) ] . [ r ( u , v ) ] 2.6.2 

Transposing both sides 

T T T -IT 

[J(v)] -fr] = [ £ (u , v) 3 .[[J(u)] ] 2.6.3 

T -1 

premultiplying both sides by [fJ(v)] ] 

X X “1 

[r] = C C JCv) ] ] . [ r(u, v) 3 -C[J(u)]'] T 2.6.4 

or 

[r] = [J(u)i“ 1 .(r(u,v)].[[tJ(v)] /r ]" 1 ] T 2.6.5 

where [r] is the required matrix of Bezier control points of 
the surface defined by [£(u,v)]. 

2 . 7 B-Spline control point* of a surface . 

A procedure similar to that described in section (2.6) 
for Bezier surfaces may be used to obtain the b-spline 
control points of a given surface. 
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A b-spline surface is defined by equation (2.4.7). 
Rewriting this equation in the matrix form 

[N(u) ] . [r] . [N(v)] = [ r ( u , v) ] 2.7.1 

following the same procedure as used in section (2.6) 

[r] = [N(u)f 1 .[r(u,v)].[[[N(v)] T ]" 1 ] T 2.7.2 

where [r] is the required matrix of b-spllne control 
vertices . 

2 . 8 Intersections of surfaces with arbitrary planes , 

A heuristic approach has been adopted to get the 
intersection curve between a surface and a plane. The 
overall procedure has been outlined below : 

(1) discretization of the surface into rectangular regions 

(2) setting up of transformation matrix such that the 
cutting plane becomes the z=0 plane 

(3) transformation of surface data points using 
transformation obtained in (2) 

(4) identification of dicretized regions through which the 
cutting plane passes 

(5) calculation of Intersection points in each of the 
regions identified in (4) 

(6) setting up of inverse of transformation obtained in (2) 

(7) drawing of cutting plane and intersection curve. 
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2 - 8 . 1 Discretization o f the eurfa c e . 


The surface points are taken in order of increasing. u 
and v and are discretized into rectangular regions (not 
necessarily planar) in the u-v plane. The number of regions 
is a function of the complexity of the surface. It is 
assumed that the surface does not contain any 
interpenetrating areas . 


This is a preliminary step both for intersections as 
veil as for hidden surface display. 


2 . 8.2 Transformation . 

To ease up the calculations required for intersection 
calculations a transformation of the coordinate system to 
the cutting plane is required. This consists of translating 
the object coordinate system to a point on the cutting plane 
and then aligning the z-axia with the normal to the cutting 
plane. This results in searching for those regions of the 
surface where z=0 lies. The steps taken in achieving this 
transformation is given below ; 

Let the three non-col 1 inear points defining the cutting 
plane be and The normal to the plane from the 
origin of the coordinate system is then given by 
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2 . 8 . 2 . 1 


■ - * Sa ~ ^ 

I ( E3 - £l) x (£j - El ) | 

the distance of this plane from the origin is given by 


d ■ Q • El 2 . 8 . 2. 2 

the point at which the normal intersects the plane is 

* ■ d Q 2. 8. 2. 3 

Setting up a translation transformation so that the origin 
of the coordinate system gets mapped on to the point defined 
by the vector r , the translation matrix is 


T 


l 

0 

0 


-r 


x 


~ r y -r z 


o 

0 

0 

1 , 


2 . 8 . 2. 4 


Since the normal to the cutting plane is to be aligned to 
the z-axis of the coordinate system, the z-axis direction is 
defined by n. Taking an arbitrary vector lying on the 
cutting' plane, say b = r) and assigning this to the y- 

axis direction, and the vector product a of b and the normal 
vector, as the x-axis direction, a rotation matrix R can be 
s et up 

Bw by ^y ^ 

r - “y b y n y ? 

i "«? ; i 

2 . 8 . 2. 5 


27 



Post multi plying equation (2. 8. 2. 4) with (2. 8. 2. 5) results in 
the required transformation matrix T 

a X b x n X ® ' 

Qy by ny 0 

a 2 b z n 2 o 2 . 8 . 2 . 6 

~ r x “ r y ~ r z 1 < 

The inverse of this matrix is also required for later 
back-transformation to the original coordinate system. 



2.8.3 Calculation of Intersection points . 

Setting up a parametric interpolation function, let the 
z-values of a rectangular region of a surface saddle a zero, 
with f i z and f^ z the z-values at the two vertices between 
which there is a intersection point, then 


(l-t).f lz * t.f 2z - 0 2. 8. 3.1 

where t is the parametric variable lying in the interval 
0 < t < 1 2. 8. 3. 2 

therefore 


t 



f 2z _f lz 


2 . 8 . 3. 3 


Then using the value of t obtained in the above equation the 
x and y-values of the intersection points are given by 
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X - (l-t).f lx ♦ t.f 2x , and 

y - ♦ t.f 2y 


2 . 8 . 3 . 4 
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2 . 9 Hidden Surface Diaplay . 


The rectangular vertices obtained after discretization 
of the surface are in the object coordinate sytem, and 
therefore they have to be transformed to the viewing plane 
coordinates before any hidden surface algorithm can be 
appl i ed . 

The surface data points are first transformed to the 
viewing plane. It can be assumed that the rectangular 
regions are almost planar and that they do not 
interpenetrate. The average depth of each rectangle is 
therefore the mean of the z-values of the four vertices. 

Using a bubble sort the rectangles are arranged in 
depth-order with the rectangle farthest from the viewer 
appearing first in the list of sorted rectangles. 

The sorted list of rectangles is then accessed and 
and each is filled with a colour selected by the user. The 
rectangle drawn last covers up the already drawn rectangles, 
thereby creating a hidden surface view. Uith this technique 
the hidden line problem is trivially solved if ‘the fill- 
colour used is the background colour. 


CHAPTER III 


SURF-3D 


3 . 1 Overview of the system, 

SURF-3D ia a graphic interactive, menu-driven software 
for the design and display of three-dimensional surfaces. 
The formulations described in sections 2.2, 2.3 and 2.4 have 
been imp-1 ement ed in this package, thereby facilitating the 
design o-f spline-blended, Bezier and B-spline surfaces. A 
comprehensive set of viewing parameters have been introduced 
to enable the user to visualize the surface in any 
orientation, either in axonometric or perspective 
projections. In addition hidden surface and hidden line 
views can also be obtained to gain more realism. 

The surfaces can be modified by changing the positions 
of control points which may be graphically ’picked’ from the 
screen for input to the program. 

Facilities for general curve generation for later use 
in designing swept, axl-symmetr ic surfaces have also been 
included as a special feature of the software. 

Input data may either be created on-line or may be 
given as input from a file, thereby making it easy for 
available surface data to be Input. 

All data regarding the control points, curve data 
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points, and surface data points can be obtained as output 
from the program in a form suitable for documentation 
purposes . 

The entire system is modular in structure and open- 
ended to the extent that additional features may be added in 
future without having to re-edit or change the package. All 
the modules have been written using FORTRAN-77 and are 
implemented in the Norsk Data 560/CX computer system. The 
program presently runs on Tektronix 4107 and 4109 graphic 
terminals. Changes may be required if the package is to be 
run on other systems. The underlying basic graphic action 
routines are built-up on the PLOT-10 Terminal Control System 
software, a product of Tektronix Inc. An entire suite of 
three-dimensional wireframe graphics routines has been 
implemented in this package. 

The data structures used are simple and are close to 
their actual mathematical equivalents e.g. a set of vectors 
has been represented as a three-dimensional array in the 
program . 

3 . 2 Int eract ive menu f eaturea . 

The system provides an interactive menu, which helps in 
going through the steps of input, design and display in a 
systematic manner. The user at any stage is given a display 
q £ all the commands available to him at that stage. The 
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main advantage of this is that the user does not have to 
commit to memory any of the commands. 

There are in all six menus viz. 

(i) MAIN 
Cii) VIEU 
(lii) CURVE 
Civ) FIT CURVE 

(v) SURFACES 

(vi) SU ITCHES 

Each of these above provides sets of options or 
commands which can be selected when the program displays a 
particular menu. The user, for example, may input data from 
a file using the MAIN, then set the viewing parameters using 
VIEU and then on to surface generation and display using 
SURFACES, which gives a choice of selecting among the three 
different types of surfaces presently supported, or for 
intersections and hidden surface displays. 

The selection process has been made simple. The user 
needs only to align the cursor with a particular option and 
then press the space-bar to complete the selection. The 
program gives an immediate feedback if a proper selection is 
made. 


The overall structure of the menu system is modular and 
new commands and options may be added as and when required. 
The structure of the menu system has been depicted in figure 
( 3 . 2 . 1 ) . 
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Fig. 3.2.1 Menu Structure. 


















































3 . 3 Graphics f eaturea . 


A full-fledged three dimensional wire-frame graphics 
system has been implemented as a part of the package. This 
has been built up as a super-set over the set of basic 
graphics action routines provided by the PLOTIO- Terminal 
Control System. This therefore, has resulted in efficient 
use of all the graphic facilities offered by Tektronix 4100 
series colour graphic terminals. 

The following are the features of the three-dimensional 
graphics system : 

(i) full three-dimensional windowing has been provided so 
that us er-speci f i ed coordinate values can be easily 
mapped to the screen-coordinate system. The user works 
in an environment in which he neither requires knowledge 
of computer graphics nor of computer pr ogrammning . 

(ii) three dimensional viewing transformations encompassing 
both parallel and perspective projections, including 
independent translation, scaling and rotation about 
any axis are possible 

(lii) hidden surface display usln& panel-filling of surface 
patches with user selectable colours. 


s 
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3.3.1 Dlaplay optloni, 


The three-dimensional graphics routines support full 
colour viewing transformations in the principle of a 
synthetic camera. The user can change the view parameters 
like direction of projection, center of projection or eye- 
point, view-up direction, view distance, view plane position 
and orientation etc. to arrive at a required view for 
display . 


In addition, a user may specify independent three- 
dimensional coordinate axes rotations, translations, and 
scaling transformations for viewing. 


Facilities have been provided to selectively display 
entities like points, polyhedrons or polygons, curves, 
surfaces, axes triad etc. This allows the user to specify 
explicitly what he wants to display. Points are displayed as 
identifiable marked points, with control polygons and 
polyhedrons appearing in straight lines. The surfaces and 
curves are drawn as families of curves. 


3.3.2 Viewing transformation . 

The view parameters that the user can set within the 
program are : 


(i) 


the view reference point (XR,YR,ZR); the 
translation and scaling transformations are 
about this point. (The default value is (0,0 


rotation , 
all done 
0 )) 
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cm 


the view plane normal (VPX , VPY , VPZ) is a vector which 
ia perpendicular to the view plane and pointed towards 
the view reference point. (The default is the z-axis 
( 0 , 0 , - 1 )) 

(iii) the view distance (VD) is the distance of the view 
plane from the view reference point( default 0) 

(iv) the view-up direction vector ( VPX , VPY , VPZ) ; the 

direction which should be vertical, when the viewing 
transformation is applied( default being the y-axis 
( 0 , 1 , 0 )) 

These parameters allow the user to select a particular 
view of the surface. 

For parallel projection, the user must specify the 
direction of projection ( PX , PY , PZ )( de fault (0,0,1)), whereas 
in the case of perspective projections, the user has to 
specify the center of projection point or eye position 
(XC, YC,ZC) . 

Once all of these have been set by the user, the 
following steps outline the transformation procedure : 

Step 1 .. Translation of the object coordinate system to 
the view plane origin; 

Step 2. Rotation about the x-axis to place the view plane 
normal in the x-z plane; 

St ep3 . Rotation about y-axis to align the z-axis of the 
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coordinate system with the view plane normal; 

Step 4. Rotation about z-axis to align the view-up vector 


with the y-axis of the coordinate system; 

Step 5. If perspective projection is specified, then 
transformation of the eye point to the view plane 
coordinat e system; 
else 

if parallel projection, then transformation of the 
parallel projection vector to the new coordinate 
system . 

3 . 4 Surface modules 

The following surface modules have been implemented: 
(i) Spline surface 
(ii - ) Bezier surface 
(iii) B-spline surface 


3.4.1 Spline surface . 

Input to this module consistes of 

(i) the rectangular array of control points 

(ii) number of intermediate points per span 

Assuming free ends for each curve on the surface unless 
otherwise specified, a general spline curve segment between 
points r | and r^ +i and using equation (2.2.3) 

# 

rj(w) - rj.(l-3w 2 +2w 3 ) + rj.,. j.(3w 2 -2w 3 ) ♦ r* j.(w-2w 2 +w 3 ) 3.4.1 

♦ I't+iC-w^w 3 ). 
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Differentiating with respect to v 

l'l(w) - £^-6^+6^) f rj + jCGw-Sw 2 ) +£' j (1-4W+3W 2 ) 

+ l'ii-lt"2w+3w 2 ) 3.4.2 

and again 

I|(w) “ £j(-6+12w) ♦ Ij + ^(6-12w) ♦ r’j(-4+6w) 

+ r’ 1 + 1 (-2^6w) 3.4.3 

the slopes rjs are unknown; substituting v=0 and v=l in the 
equation (3.4.3) 

rffO) - -&r t + 6i 1 + 1 - 4r*j - 2£‘ 1 + 1 3.4.4 

£|C1] - 6£j - 6£ j + 1 ♦ 2r‘ j ♦ 4£' i + 1 3.4.5 

Since second derivative continuity is assumed 

r",Cl) - r 1 + 1 (0) 3 ' 4 - 6 

therefore from equations (3.4.4) and (3.4.5) 

r' t ♦ 4r’, + 1 ♦ L\+2 “ _3£ i + 3£ l +2 3.4.7 
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for free ends setting i=l 


2r'i + r '2 = -3r j + 3£2 


3.4.8 


similarly for the other end of the spline setting i=n _ l 
I'n-l + 2 l‘n * - 3 En-l + 3r -n 

3.4.9 

Using equations (3.4.8) and (3.4.9) for the ends and 
equation (3.4.7) for the other points of the spline 

[A].[r’] = [ B ] . [ r ] 3.4.10 

where [A] is n x n matrix and [B] is a matrix with n rows 
Therefore 

-1 

[£’ ] = [A] . [B] . [r] 3.4.11 

which gives the required slopes in one direction at each of 
the points on the curve. 

Substitution of the points and slopes in the equation 
(3.4.1) and stepping through unit intervals of the 
parametric variable w for each span, a curve on the surface 
is generated. 
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3 . 4.2 Bezi er surface . 


Input to this module consists of : 

(i) the rectangular array of control points 
(11) number of sections on the surf ace(nsec) 

The maximum size of the rectangular array of control 
points allowed by the program is 20 x 20. The number of 

sections (nsec) determines the number of surface points that 

are to be calculated. The number of surface points is given 
by (nsec+l) x (nsec+l). The maximum number of surface points 
that can be calculated by this module is 50 x 50. 

The following steps outline the algorithm for 

generating a Bezier surface. 

Step 1. Set number of points on each curve (np) = naec+l 
Step 2. Initialize u and v parametric values required for 
generating np x np points on the surface: 
ustep = 1 /(np-1) 
vstep = l/(np-l) 
for i=l to np do 

u[i] = ( i-1 ) * ustep 
v[i] = (i-1) * vstep 

od ; 

Step 3. Repeat steps 4 to 6 for each value of u. 

Step 4. Find blending function values using equations 
(2.3.1) and (2.3.2) 

Step 5. Calculate intermediate curve-defining points 



using equation (2.3.3) 


Step 6. Repeat steps 7 to 8 for each value of v 
Step 7. Find blending function values in v-direction 
using equations (2.3.1) and (2.3.2) 

Step 8. Calculate surface points using equation (2.3.4). 
Step 9. Draw surface. 
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3.4.3 B-»pl in« surface. 


The Input to this module consists of 

(i) the rectangular array of control points 

(il) number of sections on surf acefnsec ) 

The maximum size of the rectangular array of control 
points al-lowed by the program is 20 x 20. The number of 
sections (nsec) determines the number of surface points that 
are to be calculated. The number of surface points Is given 
by (nsec+1) x (nsec+l). The maximum number of surface points 
that can be calculated by this module is 50 x 50. 

The following steps outline the algorithm for 
generating a B-spline surface. 

Step 1. Set number of points on each curve (np) = nsec+1 
Step 2. Calculate the knot values in both u and v 
directions using equations 2.4.4 and 2.4.5 
Step 2. Initialize u and v parametric values required for 
generating np x np points on the surface: 
ustep = (t*^ )^/(np-l) 
vstep = ( t^k ) v / (np-1 ) 
for 1=1 to np do 

u [ i ] = (i-1) * ustep 

v[i] = (i-1) * vstep 

od ; 

Step 3. Repeat steps 4 to 6 for each value of u. 

Step 4. Find blending function values using equations 
(2.4.1) and (2.4.2) 
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Step 5. Calculate intermediate curve-defining points 
using equation (2.4.3) 

Step 6. Repeat steps 7 to 8 for each value of v 

Step 7. Find blending function values in v-direction 

using equations (2.4.1) and (2.4.2) 

Step 8. Calculate surface points using equation (2,4.7). 
Step 9. Draw surface. 
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CHAPTER IV 


CASE STUDIES 

4 . l Examples o f spline blended surfaces . 

The surface passes through all the given control points 
in a spline-blended surface. Therefore a change in the 
position of a control point will lead to change in the 
entire surface because second derivative continuity has to 
be maintained throughout any curve on the surface. The 
alterations that occur in a spline blended surface are 
usually unpredictable and therefore not suitable for 
modifications- It is advisable only to input available 
surface data using this representation and do all 
modifications after a conversion either to a Bezier or a 
b_spl ine surface. See figure 4.1. which is a graphic output 

of the design package. 
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Exampl eg of Bezier surfaces . 


. 2 


Bezier surfaces only approximate the control points, as 
stated before. Only the four corner control points lie on 
the surface. The blending functions used for Bezier surfaces 
of different orders give it the convex-hull property. The 
surface is contained entirely within the polyhedron of 
control points. Modifications to positions of control points 
result in the surface getting altered more near the control 
point and in varying degrees over the rest of the surface. 
However, the changes are quite predictable, and with a 
little practice a user can achieve the design he requires. 

The main drawback of a Bezier surface is that the order 
of the surface is explicitly controlled by the number of 
control points, and therefore the degree of the resulting 
surface is usually unalterable. Another shortcoming of this 
scheme is the global nature of change that occurs when a 
control point is modified. This does not allow a designer to 
make fine adjustments to the surface. See figure 4.2. 
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4 . 3 Examples of B-apline surfaces . 


A b-spline surface behaves in the same way as a Bezier 
surface, only it follows the control polyhedron more 
closely. The level of continuity of the interpolated surface 
can be explicitly controlled by the designer and therefore 
its complexity. Also the modifications are of local nature. 
These two properties makes the b-spline representation one 
of the most popular schemes for surface design. See figure 
4 . 3 . 
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4 . 4 Examples of Intersections . 


Figure 4.4.1 and 4.4.2 show intersection curves 
resulting out of a surfaces being cut by a plane. The 
Intersection curve may be obtained for a surface oriented in 
any view. 




A. 104 18.6 


51 


Figure 4.4.1. Intersection 



Figure 4.4.2. Intersection 
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4 . 5 Examp 1 es o f editing surfaces . 


Figure 4.5.1 shows a spline surface. The position of 
one of the control points was changed and a surface 
refitted. Note the alteration in the overall surface in 
f igur e 4.5.2. 

Figure 4.5.3 is a Bezier surface, which when subjected 
to editing resulted in the surface in figure 4.5.4. 

The B-spline surface depicted in figure 4.5.5 was edited 
in a similar fashion. Note the local change in the surface 
.after editing in figure 4.5.6. 
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L r\ 


F igur e 4.6.1, 


4.6.2 and 4.6.3 are illustrations of 


swept axi-symmet r ic surfaces, which can be designed by the 
package . 
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Figure 4.6.3. An axi-symmetric component (360 deg) 




4 . 7 Examples of hidden surface removal . 

Figure 4.7.1 is a program output of a surface with 
hidden faces removed by the algorithm described in chapter 
III. Figure 4.7.2 is another example. The same algorithm can 
be used to generate a hidden line view of a surface. 






CHAPTER V 


CONCLUSIONS 

5 . 1 Technical Summary 

Designers and draftsmen use the classical 'lofting' 

techniques to design surfaces. These methods are based on 

cross-sectional details obtained from templates. The entire 

procedure of design, keeping aesthetics and/or aerodynamics 

* 

in view, is time consuming and is of a specialized nature. 
Proficiency in this craft takes years of experience in 
surface design and engineering and the methods have very 
limited applications in the field of free-form sculptured 
surfaces. 

The emerging technology of computer graphics coupled 
with the advances made in computational and differential 
geometry, offers a powerful tool to designers to enable them 
to solve many of their problems in design and display of 
surfaces. In the present work, an attempt has been made to 
bridge the gap between mathematical formulations of surfaces 
and engineering design. 

A package, SURF-3D , for design and display of three- 
dimensional surfaces has been developed and implemented. 
Design of surfaces, in this software package, is based on 
the spline, Bezier and b-spline formulations, atreatmeht of 
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which is given in chapter II. These schemes not only allow 
a designer to input and modify available surfaces, they also 
accommodate design and display of surfaces in a graphic 
interactive environment. Special routines for visualizing 
the design in parallel or perspective projections has been 
developed. 

A number of design problems has been tried out. 
Illustrations of some of the case studies are given in 
chapter IV. 

The salient features of the anlytical as well as the 
computational work are as follows: 

(i) A surface generation module has been developed, which 
can cater to spline, Bezier, and b-spline schemes of 
representation. The algorithms developed are 
computationally fast and therefore, are useful in an 
interactive environment . 

(ii) A editing module for creation and modification of 
control points has been implemented. The routines use 
graphical techniques to specify points. 

(iii) A view module has been devised to cater to three- 
dimensional transformations for both parallel and 
perspective projections. To ease the burden of th€ 
user in specifying viewing transformations, the 
routines simulate a synthetic camera. In addition, 
independent viewing transformations like rotations 
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translations, and scaling have been provided. 

A simple hidden surface removal algorithm has been 
developed. It is based on z-depth sort. It is fast 
enough for interactive work. 

(iv) A curve module has been Installed which is used for 
generating swept, axy-symmetr i c surfaces. Both, 
editing and the view modules can be used on these 
sur faces. 

(v) A conversion module has been implemented for inter- 
conversion between surface types. This is primarily 
meant to ease modifications. 


5.2 


for further work. 


(1) The present work deals with single instances of spline, 
Bezier and b-spline surfaces. In practice, blending 
between surfaces of different types are required. It Is 
necessary to take Into account such applications. 


C 1 1 ") Besides intersections of surfaces with planes, one 
often needs to get intersections curves of two 
surfaces. Nod i f I cat i ons and/or additions to the current 
intersection routines, or development of a general 
purpose algorithm for both planes and surfaces may be 
requlr ed . 


(iii)A very wide scope of development, exists in the areas of 
boolean set operations like cutting, joining. 
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di f f erence 


and 


intersections of 


surface(s) . 


In 


addition, interfaces to other solid modelling systems 
may be thought of. 

Civ) Standard exchange specifications like IGES and MAP, for 
transferring geometric data, can be implemented to 
enhance the program capabilities to pass as well as 
take as input, data to and from other design systems. 

(v) For manufacturing applications, numerically controlled 
machine interfaces may be developed. The design modules 
may be enlarged to take care of constraints of die- 
sinking and die-design applications. 

(vi) A more versatile hidden surface removal can be thought 
of. The current algorithm does not work for complex 
surfaces, especially those with interpenetrating areas. 
In addition a full-colour shaded image display 
algorithm may be developed to give more realism to the 
displays provided by the package. 
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